Non-equilibrium phase transitions in tubulation by molecular motors 
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The extraction of membrane tubes by molecular motors is known to play an important role for the 
transport properties of eukaryotic cells. By studying a generic class of models for the tube extraction, 
we discover a rich phase diagram. In particular we show that the density of motors along the tube can 
exhibit shocks, inverse shocks and plateaux, depending on parameters which could in principle be probed 
experimentally. In addition the phase diagram exhibits interesting reentrant behavior. 
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Molecular motors play a fundamental role in intracellu- 
lar traffic Q, being responsible for the transport of vesi- 
cles and the extraction of membrane nanotubes (T\. The 
latter phenomenon is of particular interest as it requires a 
cooperative effort between many motors. This remarkable 
collective behaviour has been demonstrated in vitro only 
very recently [2|, triggering much interest in the features 
of the "tubulation", with a particular focus on the dynam- 
ics of the tip region. Objects of study include the conditions 
for the formation of tubes, their velocities, the load exerted 
on and by the motors, the distribution of motors along the 
tube and the role of processivity ||2l |3j |4l |5l. In vitro, a 
tube can be created when a vesicle coated with kinesins is 
brought near a microtubule It is generally believed JSj 
that two regimes are then observed depending on the mo- 
tor density: below a critical density, the motors which bind 
to the microtubule are not able to extract a tube; above the 
critical density tubulation occurs and the motors pull a tube 
out of the vesicle, at steady velocity. During this process 
motors constantly bind and unbind from the microtubule, 
while remaining bound to the membrane (see figure [TJ. In 
this regime, the density of motors is predicted to be flat, 
with some structure near the tip region O. The critical 
density and the velocity of the tube have been shown to be 
very sensitive to details of the tip region, such as the num- 
ber of motors clustered there and their coordination Q. 

In this work, on the other hand, we focus on the regime 
where tubulation is established and study the collective be- 
haviour of motors in the bulk of the system. We consider a 
generic model of two coupled lattices representing bound 
and unbound motors. By accounting for excluded volume 
interactions, neglected in [3|, we discover a richer phase 
diagram than previously expected. The tubulation regime 
now divides into two different phases with re-entrant tran- 
sitions between them. Both phases could in principle be 
accessible to experiment by control of the vesicle density. 
As we show, the phase diagram is governed by the bulk dy- 
namics and the effective tip velocity; it is thus insensitive 
to the precise details of the dynamics of the tip region. 

We first describe the phenomenology predicted by our 
study. The motor density profile comprises two plateaux 
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Figure 1: Illustration of tube extraction. Molecular motors are 
attached to the membrane and can bind and unbind from the mi- 
crotubule. Two frames of reference are used in the text: the 'tip 
frame' is comoving with the tip, with site labels increasing to- 
ward the vesicle; the 'lab frame' where the vesicle is stationary 
and the site labels increase towards the tip. 



emerging from the tip of the tube and the vesicle, respec- 
tively. These plateaux meet in the bulk of the system which 
leads to a discontinuity - a kink - in the density profile. The 
system can be in two different phases, illustrated in figure 
[2] (/) a kink phase, in which the tip density is either larger 
or smaller than that of the vesicle, the two plateaux being 
connected accordingly by a shock or an inverse shock in 
the bulk of the system, and (ii) a tip phase where the kink 
travels toward the vesicle and localizes in its viscinity, thus 
yielding a constant density profile corresponding to the tip 
density. Apart from a carefully chosen set of parameters, 
the kink is never at rest, always travelling away from the 
tip of the tube and either toward (tip phase) or away from 
the vesicle (kink phase). In the latter case, the kink moves 
away from both boundaries, which is possible because the 
tube is extending. The transition between the two differ- 
ent phases, and also the transition from 'shock' profile to 
'inverse-shock' profiles within the kink phase, can be trig- 
gered, for instance, by changing the value of the vesicle's 
density (see the phase diagrams in figure [2|. Furthermore, 
the phase diagram is re-entrant: by continuously increasing 
the vesicle's density, one can go from the kink phase to the 
tip phase and back again into the kink phase. 

In non-equilibrium statistical physics, shocks play an im- 
portant role for driven lattice gas models |i7j| but the phe- 
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Figure 2: Possible phase diagrams depending on the value of 
the ratio d/a of detachment to attachment rate. The axes rep- 
resent the densities (rt, tv) of bound motors in the tip and vesicle 
plateaux, (a) 1/4 < d/a. The system presents either shock or in- 
verse shock profiles, (b) 1/8 < d/a< 1/4. The tip phase appears 
inside the shock region, (c) d/a < 1/8. The boundary of the tip 
phase moves into the inverse shock region. In cases (b) and (c), 
reentrant transitions are possible, e.g. along the blue line. 
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Figure 3: Coupled lattice model for the bulk dynamics. The sites 
are labelled in the lab frame. The bottom and top lattices repre- 
sent bound and unbound motors respectively. Possible transitions 
of motors in the bulk are illustrated by arrows in the left panel. 
The right panel shows how extension and retraction of the tube 
drag the unbound motors in the bulk. 



nomenology described here differs from previously studied 
biophysical traffic problems [9|; reentrance is unusual and 
previously observed inverse shocks required slow particles, 
static defects or special current-density relations [i8 |. 

Definition of the model. The microscopic details of the 
model are as follows (see figure [3]l. We consider two cou- 
pled one-dimensional lattices, for bound and unbound mo- 
tors, which extend from the vesicle to the tip of the tube. A 
motor bound to the microtubule steps toward the tip of the 
tube at rate p, provided the arrival site is empty. Each site of 
the unbound motor lattice accounts for a whole perimeter 
of the tube, which in experiments exceeds 100 nm [2 | and 
can contain many motors. For sake of clarity we thus ne- 
glect the exclusion on this lattice and assume that unbound 
motors diffuse freely at rate D. Partial exclusion could be 
taken into account: for realistic values of the parameters 
(see below), it does not modify qualitatively the phase dia- 
gram and just obscures the algebra [ 10 1. Finally, motors at- 
tach at rate a to an empty site and detach at rate d from the 
microtubule. A complete description of the tube dynam- 
ics would also include the details of the dynamics in the 
viscinities of the tip and the vesicle. As shown numerically 
in lO, these details are important for establishing the con- 
ditions for tubulation. However, as demonstrated below. 



the form of the phase diagram is insensitive to these details 
and relies on the fact that the tube has a well defined mean 
velocity during tubulation. We thus posit that extension 
and retraction events occur with rates w+ and , yielding 
a tube velocity v^p = —v~ . Also, since the viscosity of 
the membrane is two orders of magnitude larger than that 
of the buffer [3|, the unbound motors are dragged every 
time the tube extends or retracts (see figure [3]l. Under these 
conditions our results encompass a whole class of models 
for the tube dynamics including, for instance, those con- 
sidered in [6|. We now derive the different phases within a 
mean-field analysis. 

Mean-Field (MF) Theory. In the analysis that follows it 
will be useful to consider two distinct frames of reference: 
the lab frame, where the vesicle is stationary and the site 
labelling starts at the vesicle and increases toward the tip 
of the tube; the tip frame, which is co-moving with the tip 
of the tube, and where the site labelling starts at the tip and 
increases toward the vesicle (see figure[T]). In the tip frame, 
the mean-field equations read 

= -J^ + Jl,+K, ■ a, = -J^ + Jl^-K, . (1) 

Here, and a.t are the average occupancies of bound 
and unbound motors at site i. The current of bound mo- 
tors moving between sites i + 1 and i is given by = 
—p{l—Ti)Ti^i+v'^Ti—v'Ti^i whereas the current of un- 
bound motors reads J" = — Z)[(Ti+i — 0"^]. Both currents 
are defined to be positive in the direction of increasing i, 
i.e. when they transport motors away from the tip. Also, 
in the tip frame, extension and retraction of the tube affects 
the bound motors throughout the lattice, whence the con- 
tribution of v+ andv^ to J^. Last, Ki = aai{l — Ti) — dvi 
is the flux of motors between the two lattices. The counter- 
part of these equations in the lab frame (i increasing toward 
the tip) is easily obtained: 

f. = -j':+jU+K, ; a, = -j: + jt,-K,, (2) 

where = -D{a,+i - CTj) + v+ai - v^a,+i and j\ = 
pTi{\ — Ti+i) are the unbound and bound currents. Here 
also, currents are positive in the direction of increasing i. 

As we now show, the velocity of the tip of the tube 
selects plateau densities of bound and unbound motors, 
which we call Tt and at respectively. At the other end, 
the density of motors on the vesicle selects in general dif- 
ferent plateau densities which we call Ty and a^. To derive 
the steady-state plateau densities, we assume constant Tt^v 
and (Tt,y in either ([T]l or Q. This yields a zero flux between 
the two lattices Ki = 0, implying for any pair of plateau 
densities r and a 

a = dT/[a{l - r)] . (3) 

Let us first consider the tip plateau values Tt and at us- 
ing Eq. ([T). Adding upper and lower lattice contributions 
yields a conservation equation for the total flux 

Ft = + = . (4) 
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The total flux Ft flowing through the tip plateau in the tip 
frame has to equal zero as nothing can move to the left of 
site 1 . Using the explicit expressions of and J" in Q 
and relation ([3]) yields 

, vup d f p 

Tt = l ; CTt = - — 

P a \Vtip 

Note that we do not specify equation (or dynamics) near 
the tip region. Solving such equations will give (generally 
complicated) relations between the rates at tip region and 
values of v^p, Tt and at while leaving Eq. |5] unmodified. 
They will therefore not influence our results. 

The vesicle plateau, however, is determined by the den- 
sity of motors on the vesicle and the details of the nearby 
dynamics. While such equations can be solved for specific 
models, to analyze the phase diagram it is enough to know 
that Ti, can take any value between and 1. 

Tip and vesicle plateau densities are typically different 
which suggests the possibility of a kink phase with shock 
and inverse shock profiles when Tt > r„ and Tt < Ty 
respectively. Generally, the kink is not at rest and this 
phase disappears if it propagates to either end of the sys- 
tem. To analyse this, we consider the kink velocity in the 
tip frame, , and in the lab frame, vf. Conservation of 
mass implies = {Ft - F^)/{pt - p^), where Ft.y 
and pt^v = Tt.v + <7t,t) are the total fluxes and densities 
to the left and right of the kink, in the tip frame. F^ reads 
Fy = —p{l—Ty)Ty+VtipTy. Using Q and ([5]l to eliminate 
(Tf_„ and Vtip, we obtain Vt in terms of „: 

pTy{\ - Tt){l - T^ 



{l-Tt){l-Ty)+d/a ■ 



In the lab frame, the kink velocity is 
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pTy{l - Tt){l - Ty) 

{l-Tt){l-Ty)+d/a • 
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Since are smaller than 1, Vt is necessarily posi- 
tive, i.e. the kink always propagates away from the tip. 
Transposed in the lab frame, this means that the kink never 
catches up with the tip. However, f f can be negative, i.e. 
the kink may not propagate away from the vesicle. The 
tip phase indeed occurs when the kink is localized at the 
vesicle and the density then equals that of the tip plateau, 
except in a boundary layer close to the vesicle. 

The tip phase thus requires < 0, which reads (1 — 
Ty){l — Ty — Tt) + d/ a < and can only be satisfied if 

Tt > 2^/d/a . (8) 
The system is then in the tip phase for Ty G [r_ ; r+] , where 



r± = 1 
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When Ty is not in [t_ ; r+J or condition ([8]l is not met, the 
system is in the kink phase, presenting shock when Tt > Ty 
or inverse shock when Tt < Ty. Note that for G [0, 1], 



one always has r_ > and r+ < 1. The phase diagram 
thus always exhibits reentrance if d/a < 1/4, i.e. there 
always exist values of Tt for which a continuous increase 
of Ty drives the system from the kink phase into the tip 
phase and back into the kink phase. We present the various 
possible phase diagrams in figure |2] 

Numerics. In order to validate the theoretical predic- 
tions, we now turn to the simulation of a concrete model 
within the class considered here. The bulk dynamics has 
already been described (see figure |3]l and we now specify 
dynamics in the viscinity of the tip and the vesicle. Our in- 
terest lies in verifying the phase diagram and to this end we 
choose a particularly simple model. The vesicle is repre- 
sented by reservoirs of bound and unbound motors of den- 
sities Tg and (Jq. For simplicity, they are chosen to satisfy 
([3]) so that there is no flux between them ifTTI . At the other 
end, the tube extends at rate 7 by one lattice site if a bound 
motor occupies the site next to the tip. When this happens, 
the new unbound site next to the vesicle is equilibrated with 
the reservoir of density aQ. If a bound motor does not oc- 
cupy the site next to the tip the tube retracts with rate p. At 
this site motors can still attach and detach with rate a and d 
and unbound motors can hop toward the vesicle, with rate 
D. At long time, this tip dynamics yields average extension 
and retraction rates = and = p{l — Ti). 

Models aiming to predict the threshold for tubulation 
and tube velocity would require modified attachment and 
detachment rates at the site closest to the tip and should 
account for backward stepping there. However, as noted 
above, the phase diagram depends only on the tip velocity 
and not on further details of the tip dynamics. 

We now consider the results of continuous time simula- 
tions of the model. In figure |4] a typical shock profile and 
its dynamics are presented. Note the quantitative agree- 
ment with the predicitions of the mean-field theory. 

We have also verified the general structure of the phase 
diagram and present results for the blue line indicated in 
figure 2 The different profiles observed are presented in 
figure 5 From previous lattice gas studies, one would ex- 
pect the inverse shock to smooth out through a rarefaction 
fan. Here, on the other hand, we checked numerically that 
its relative width vanishes in the large time limit; it is thus 
stabilized by the interaction of the two lattices |10|. The 
figure also shows that the tip velocity is independent of Ty, 
once tubulation is established. A mean-field analysis |[TOl 
suggests that this holds for generic local tip dynamics. 

To quantify the transition, we define A = ^ . Ti/L — 
Tt, where L is the length of the tube. This compares the 
average mass of bound motors in the system with that of a 
putative tip phase. The parameter A is non-zero in the kink 
phase and zero in the tip phase. An example of re-entrance 
is shown in figure [5] Starting with Tq close to 0, we see 
that A is negative, vanishing at Tq = r_ where the system 
enters the tip phase. Further increase of Tq above Tq = r_|_ 
drives the system back into the kink phase, in the inverse 
shock region, and A becomes positive. 
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Figure 4: Shock profiles obtained from numerics (p — 1, d — .08, 
a = I, D — I, fi = 9, — 0.45, tq — 0.2). Data are averaged 
over a short time window (At = 1, 000) and then over 100 simu- 
lations. Left panel: Bound (red) and unbound (blue) motor den- 
sities along the tube, in the tip frame, at time t = 60, 000. Black 
lines correspond to MF predictions ([5]l, with vup obtained from 
the simulation. Right panel: Density of motors (black) at fixed 
time intervals (t = 4.10'*; S.IO'*; 12. lO'^) in the lab frame. Tubes 
are color-coded, with light and dark blue representing low and 
high density. Note that the tip moves faster than the shock. 
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Figure 5: Left-panel: Different profiles are illustrated as in fig- 
ure|4) but with t = 120,000, d = 0.11 and tq = 0.3, 0.6, 0.9. 
Note that the length of the lattice, hence the tube velocity, is in- 
dependent of TQ. Right-panel: Numerical value of A (crosses) 
along the blue line of figure [2] When tq goes from to 1 , the 
average density of bound motors is first lower than that of the tip 
plateau (A < 0, shock), then equals it (A = 0, tip phase), and last 
overcomes it (A — 0, inverse shock). The blue lines correspond 
to the predicted boundaries of the tip phase (|9]l. There is thus 
a reentrant phase transition from the kink phase to the tip phase 
and back into the kink phase, whose boundaries are accurately 
predicted by the MF theory. 



Conclusion In this letter we have shown that the dy- 
namics of tabulation reveals a rich phenomenology, includ- 
ing shocks, inverse shocks and re-entrant phase transitions. 
This arises from the two competing densities set by the two 
ends of the tube and the dynamics of the resulting kink de- 
termines the phase structure. This picture is substantiated 
by a mean-field theory which accurately predicts the phase 
diagram, as checked by our numerics. 

Some experimental signatures of our theory are as fol- 
lows. First, the velocity of the tip should always exceed 
that of the kink. Also, once tubulation is established, the 
velocity of the tip is not sensitive to the density of mo- 
tors on the surface of the vesicle. Last, in the experiment, 
the ratio d/a ~ 0.1 is very small [6|. It should therefore 
be possible to observe the transitions to the tip phase by 
varying the density of motors on the surface of the vesi- 



cle. To explore the full phase diagram presented in figure 
|2] one needs to change microscopic rates to vary r^. Exper- 
imentally, this could be done by changing parameters such 
as the membrane surface tension or the ATP concentra- 
tion. Finally, corrections to the phase diagram due to par- 
tial exclusion among unbound motors is of order d/aNmax 
where N^ax is the maximal occupancy of the unbound lat- 
tice [Toll. Here, it would be of the order of 0.01. 
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